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Abstract. We study J/ip propagation in hot hadronic matter using a four- 
flavor chiral Lagrangian to model the dynamics and using QCD sum rules 
to model the finite size effects manifested in vertex interactions through form 
factors. Charmonium breakup due to scattering with light mesons is the primary 
impediment to continued propagation. Breakup rates introduce nontrivial 
temperature and momentum dependence into the J /if) spectral function. 



1. Introduction 

Suppression of the J/ip signal in heavy ion collisions is still considered one of the 
promising probes of our attempts to push back deconfinement and open up a four 
volume of quark-gluon plasma (QGP) from the vacuum wherein quarks and gluons 
have their confinement scale extended to macroscopic lengths. The greatest four 
volume that one can hope for in the laboratory at present is about 10 4 fm 4 since 
stable nuclei are of limited size (~ 10 fm) and system longevity is expected to be of 
the order 10 fm. Any yet, dynamical scales for strong interaction physics are at the 
femtometer and subfemtometer level, and so it is quite reasonable to approximate the 
system as being relatively large, relatively long-lived, and near thermal equilibrium. 
Then, kinetic theory of asymptotic states and vacuum cross sections is a useful tool 
to establish the baseline rates for various dynamical effects. 

The mechanism for suppression of J/ip suggested by Matsui and Satz 1 originates 
from color screening in QCD plasma which is clearly absent in hot hadronic matter. 
Screening would effectively "partonize" the bound state mesons (the QCD analog of 
ionization in atomic physics) into open subsystems with colour (quarks, antiquarks and 
glue). For charm, this possibility brings forward a probe of open (D mesons) versus 
hidden (J/ip). Open charm decays mostly semileptonically, and will therefore be 
absent in an invariant mass window of say, muon pairs near the J/ip mass. Therefore, 
a partial or total absence of a strong J/ip peak could indicate the presence of QGP. 

In order for the probe to be useful, hidden charm must be relatively unreactive to 
hadrons to survive the hadronic stage of the fireball. For instance, we do not expect 
suppression of the (/> to be a QGP signal because the (p interacts with light mesons such 
that its mean free path is of the order 2-5 fm[21El- We ask, could the same effects 
spoil charm- J / tp that most likely spoil strangeness-^ from utility? Our purpose here 
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is to assess within an effective field theory the interaction cross sections of J /ip with 
light mesons (pions, rho mesons, etc.) and to further estimate the breakup rates. The 
hope is that such activity might shed some light on the extent to which suppression 
of the J/ip as a plasma indicator is emulated by purely hadronic effects. 

2. Chiral Lagrangian 

Quantum chromodynamics at low energies is particularly challenging. Lattice 
technology, for example, has not yet advanced to the stage where one can do dynamical 
studies at finite temperature. We instead resort to very successful methods of effective 
Lagrangian descriptions for the hadronic degrees of freedom. While the symmetry 
arguments typically introduced for two light quark flavours are not as clearly applicable 
when generalizing to three and four flavours, we indicate, as is usually done, that 
symmetry breaking effects are relegated entirely to physical mass eigenstates, leaving 
the interaction completely symmetric from the beginning. The guiding principles are 
therefore chiral symmetry, unitarity, and gauge invariance. 
We start with a nonlinear sigma model 

C= -^fd^Ud^, (1) 

where U — exp(2i<f>/F v ), and <p is the pseudoscalar meson matrix appropriately 
generalized to include charm, thereby becoming 4x4, and — 135 MeV is the pion 
decay constant. A chiral covariant derivative d^U —> D^U introduces vector and axial 
vector fields and their respective coupling to pseudoscalars. Kinetic energy terms for 
spin 1 fields and generalized mass terms of the left- and right-handed fields and 
A^ arc introduced. Then, the axial vector fields are gauged away leaving to lowest 
order a set of interactions for light plus heavy pseudoscalar and vector mesons. The 
interactions are compactly written as0] 

2 

£ int = igTr (p M [d^, 4>]) ~ yTr ([0, p"] 2 ) +igTr {d»p u [/A fT]) 

+ ^Tr(V,^] 2 ) , (2) 

where is also generalized to encompass charm and is therefore a 4 x 4 matrix of vector 
mesons. If the symmetry were perfectly realized in nature, the chiral coupling constant 
g would alone be sufficient to describe the interaction strengths universally. We relax 
that requirement here and appeal to experimental results and model calculations via 
QCD sum rules for fitting the coupling constants at each three-point vertex. Gauge 
invariance then restricts the four-point functions. 

A recent measurement of T(D* ± — > Dir) — 96 ±4 ± 22 keVjKj fixes the value 
of the coupling constant go'D-K = 6.39. Coupling of the J /if) to D D and to D* D* 
vertices is also required as input to calculate the pion-included breakup cross section 
for J/ip. Using strict vector dominance, which has been done several times in the 
past[Hl 13 03 the coupling constants would be gj/^DD = 9,j/ipD*D* — 7.7. In the 
next section we indicate an alternative approach of extracting not only the coupling 
strengths but also momentum dependence at the vertices from QCD sum rules. The p- 
induced inelastic reactions require further information via coupling constants. Vector 
meson dominance and chiral symmetry would lead to g p u d = 9 P D' d* = 5.6 6 . We 
alternatively look to QCD sum rules for information on the coupling strengths and 
momentum dependences, described below. 
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Figure 1. Lowest order Feynman diagrams for the process J/tj) + 7r — > D*D. 
As is typical in these theories, there is a i-channel (a), a u-channel (b) and a 
four-point graph (c). 



3. Form Factors 

Hadrons are clearly composites of the underlying quarks whose (effective) fields 
describe pointlike physics only when all the interacting particles are on mass-shell. 
As soon as at least one of the particles in a vertex is pushed kinematically off- 
shell, the finite size effects become important. This is of course not true for gauge 
field theories {e.g. QED, QCD) whose degrees of freedom are structureless. One 
way to account for this is to use the effective Lagrangian to construct higher order 
processes (vertex corrections, loop corrections to propagators, etc.) and to formulate 
strategics for renormalization. This has not been done up to now. However, the 
same physics has been modelled using QCD sum rules to extract the momentum 
dependence of the vertex coupling constants which then "run" with the relevant 
squared four momenta contacting the vertex [TU]. The basic scheme is to use on the 
one hand, Wilson's operator product expansion to write the vertex, and on the other, 
phenomenolgy in the form of Borel transformations to write the same function as 
linear combinations of relevant four momenta. Running coupling constants have been 
obtained as follows[TTl El EH! 



9DDJ/ip{t) 
gD'D-jr(t) 



9DDj/ii e 



[(t-16.2) 2 /228l 
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(3) 



where t is the squared four momentum of the off-shell D meson, Ad =3.5 GeV, gDDp 
= 4.4, go*DTr = 6.39, and where gDDjip = 16.4. These values are different by tens of 
percents from the chiral symmetry assumption and strict vector meson dominance. 

We next illustrate our gauge invariant method of implementing form factors for 
the channel J/ip + n — > D* + D. Notice from Fig. (J that there are three Feynman 
graphs contributing to the scattering amplitude: a ^-channel, it-channel, and a four- 
point graph. The full amplitude becomes 



M 



M?hn 



L(t) h 2 (t) + Ml h x {u) h 2 (u) + M° 4pt (-g^ AH (4) 

where the superscript "0" indicates amplitudes without form factors (without running 
couplings) modulo the coupling constants, and where the four-point amplitude is 
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evaluated only after replacing "— with I' 1 ", where 

=AgT + B (j&fi+rtti,) + C (rfrti+rttf,.) 

+ D (p%p» + p^ D ) + E (pX + P^Vd* ) ■ (5) 

Gauge invariance restricts the expansion coefficients A, . . . E, albeit not uniquely. A 
Lorentz and gauge-invariant solution is the following 

A = - h 1 (t)h 2 (t) 

B = D MjM - feiM^M 

P.J/iP ■ Ptv + Pj/4, ■ PD 

C = E = 0. (6) 

This strategy accommodates extended hadrons and is applied to all processes, 
including p-induced channels. Before coming to the results, we reiterate the features 
of our form-factor implementation and underscore the importance of using 

• form factors which are covariant, 

• unique form factors at each vertex which are appropriate for the specific 
kinematics introduced at that vertex (rather than some average momentum 
transfers over all graphs), and 

• an approach where gauge invariance is strictly respected. 

Giving up on Lorentz invariance, gauge invariance, or failing to insist on the proper 
kinematics for each form factor should be considered less reliable and therefore 
unacceptable. 

4. Results 

Cross sections for J ftp dissociation are inaccessible experimentally. It is therefore 
crucial to have reliable estimates for the magnitude and energy dependence of such 
quantities. Only then can we address the earlier question regarding the extent to 
which the J /tp is unreactive to light hadrons. We report in Fig. [21 the results of some 
of the pion and rho meson induced breakup cross sections as functions of energy. The 
near-threshold and energy dependences are nontrivial and the magnitudes of the cross 
sections are several millibarns. We note that quark-exchange models have recently 
given cross sections of the order 1 mb^31- We further note that while the magnitude 
of the cross sections here are similar to other recent meson-exchange models, the 
energy dependence is quite different owing to different form factors and different 
implementation^ tTTil ITTl IT5] . 

Reaction rates within a fireball T (a>, p) are also important indicators of spatial 
and temporal scales [21 1191 120| . Here we use the formalism from Ref.0] to write the 
J/ip spectral function as 

Aj/^(uj,p) = - 2lmDj/^(uj, p) , (7) 

where Dji^ is the scalar part of the J/ip propagator. If we neglect the difference 
between longitudinal and transverse polarizations of the J ftp in the finite temperature 
medium |21| . then we have 

p — Tflj^ — F(UJ, p) 
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Figure 2. a{J/ij> + it -> D*D + D* D) and <r(J/i/> + p -» D*D*). 



where p^ 1 = (a;, p), and i* 1 is the scalar imaginary self-energy. Then, using 

IV/v, = -1/mj/Tp Im F(p 2 = m 2 7/j/) ) , (9) 

we write 

Aj^{u,p) = )2 +m 2 T 2 > ( 10 ) 

where Tj/^, contains the vacuum width and the contributions from inelastic scattering 
with pions and rhos. 

The spectral function resulting from this calculation appears in Fig. |3| We fix the 
temperature at values of 150 and 200 MeV, which are relevant values for heavy ion 
experiments. Four different values of J/ip momenta are chosen for each. We analyze 
the spectral function this way since there is no expectation that the J ftp's will be 
thermal. The spectral function for the J/ip is clearly modified at finite temperature 
due to interactions with hadronic matter. Breakup rates fed into the spectral function 
range from 5 to 30 MeV. 



5. Conclusions 



We have applied an effective Lagrangian to describe the dynamics of J/ip with 
particular focus on breakup reactions induced by light unflavoured mesons. Off-shell 
effects at the vertices were handled with QCD sum rule estimates for the running 
coupling constants, i.e. the three-point functions. Analytic continuation to the poles 
(on-shell mesons) gave further guidance to constrain couplings. Four-point functions 
were first expanded over general Lorentz structures and then constrained with Ward 
identities. The formalism is then covariant, gauge invariant, and consistently describes 
on a graph- by-graph basis the kinematics and the off-shell vertex functions. The model 
uses effective field theory, QCD sum rules for the form factors, and Ward identities to 
constrain the overall scattering amplitudes. 

Resulting cross sections for breakup of J ftp are several millibarns, and exhibit 
nontrivial energy dependence. When imported to the spectral function, the resulting 
effects give a dramatic modification at finite temperature as compared to vacuum 
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Figure 3. Spectral function in vacuum and at finite temperature for fixed 
momenta and fixed temperatures. 



structure. This could have important consequences for studies at SPS and RHIC 
measuring dilepton signals near the J/tp mass and impact interpretation of J/ip 
suppression. 
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